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Abstract: Composite dark matter is a natural setting for implementing inelastic 
dark matter — the 0(100 keV) mass splitting arises from spin-spin interactions of 
constituent fermions. In models where the constituents are charged under an axial U{1) 
gauge symmetry that also couples to the Standard Model quarks, dark matter scatters 
inelastically off Standard Model nuclei and can explain the DAMA/LIBRA annual 
modulation signal. This article describes the early Universe cosmology of a minimal 
implementation of a composite inelastic dark matter model where the dark matter is a 
meson composed of a light and a heavy quark. The synthesis of the constituent quarks 
into dark hadrons results in several qualitatively different configurations of the resulting 
dark matter composition depending on the relative mass scales in the system. 



Contents 



1. Introduction 1 

1.1 Review of Axial CiDM Model 2 

2. Dark Matter Synthesis Prelude: Hadron Spectroscopy 4 

2.1 Coulombic Binding Energies 5 

2.2 Hadronization of Bound States 7 

2.3 Hyperfine Structure Spectroscopy 9 

2.4 Light States 11 

3. Dark Matter Synthesis and Evolution 12 

3.1 Preconfinement Synthesis 12 

3.2 Synthesis at Confinement 14 

3.3 Post confinement Synthesis 14 

3.4 Synthesis Results 16 

4. The Excited Fraction of CiDM 18 

4.1 Early Time Depopulation 18 

4.2 Late Time Up Scattering 19 

4.3 Down Scattering Discovery 20 

5. Baryon Elastic Component and Direct Detection 22 

5.1 Heavy Baryon Scattering 22 

5.2 Light Baryon Scattering 24 

6. Conclusion 25 



1. Introduction 

DAMA/LIBRA has an on-going 8.9cr annual modulation signal of their single hit rate 
[1-3]. If this signal is caused by dark matter scattering off their Nal crystals, then 
DAMA, when combined with other null direct detection searches, suggests that dark 
matter is scattering inelastically to an excited state split in energy by (9(100 keV) 
[6-9]. Inelastic dark matter (iDM) is a testable scenario with XENONIO and ZEPLIN3 
performing dedicated reanalyses of their data [4, 5] . IDM will be decisively confirmed 
or refuted in the XENONIOO, LUX, or CRESST science run this year [10, 11]. 
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Composite inelastic Dark Matter (CiDM) is a model of iDM, where the (9(100 keV) 
scale is dynamically generated by hyperfine interactions of a composite particle [12]. 
(See [13-20] for other examples of composite dark matter.) In the original CiDM 
model, the dark matter consists of a spin-0 meson, vTd, that has a single heavy con- 
stituent quark. Adjacent in mass to ttj is the spin-1 dark meson, pd? and with the mass 
scales chosen in [12], the mass splitting between and pd is 0(100 keV). The dark- 
matter origin in CiDM is non-thermal and requires a primordial asymmetry between 
the number densities of heavy quarks to heavy antiquarks. This article addresses the 
early Universe cosmology of this minimal CiDM model, calculating the abundance of 
TTd bound states relative to pd states, dark baryon states or other exotic configurations 
of the dark quarks, as well as the direct detection properties of the various components. 

The results of this article show that a wide range of final abundances of dark 
hadrons is possible, depending on the spectroscopy of the dark sector states. Some 
spectra have dark matter dominantly in the form of TTd, while other spectra have dark 
baryons dominating the abundance. In the latter case, where the dark baryons are the 
predominant dark matter constituent, the residual vTd component interacts sufficiently 
strongly to account for DAMA's signal. In all cases, exotic dark matter components 
arise in CiDM with novel elastic scattering properties - nuclear recoil events are sup- 
pressed at low-energy by a dark matter form factor. The relative abundance of pd to vTd 
is typically rip^/riT,^ ~ 10^"^. Existing searches for down-scattering (pd — )■ vTd) in direct 
detection experiments are not constraining, but may be feasible in the near future. 

The organization of this article is as follows. Sec. 1.1 briefly reviews a specific im- 
plementation of CiDM presented in [12] that will be used throughout this analysis. The 
spectroscopy of this model is discussed in Sec. 2 and the synthesis of dark meson and 
dark baryon states in the early Universe in Sec. 3. Sec. 3.4 is the primary result of this 
article and the qualitatively different results of the synthesis calculation are classified 
here. Sec. 4 addresses the upscattered fraction of dark pions and its implications for 
direct detection. Sec. 5 summarizes constraints on the dark baryon fraction that arise 
from direct dark matter searches such as CDMS and comments on the novel proper- 
ties of the elastic scattering processes. Sec. 6 concludes with the outlook and further 
possibilities. 

1.1 Review of Axial CiDM Model 

A wide variety of hidden sectors weakly coupled to the Standard Model have been 
considered in the literature, and the possibility that dark matter is charged under 
hidden sector gauge forces has received considerable recent attention [21-32]. In CiDM, 
dark matter is a bound state with constituents charged under a hidden sector gauge 
group of the form SU{Nc) x t/(l)d, where the U{l)d gauge boson kinetically mixes 
with Standard Model hypercharge gauge boson, and the SU{Nc) group condenses at 
the GeV scale. Theories in which dark matter is charged under a new GeV-scale 
gauge group, as in CiDM, predict a variety of multi-lepton signals in _B-f actor ies and 0- 
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factories [33-37], low-energy upcoming fixed-target experiments [38-41], and distinctive 
astropliysical signatures [42-44]. 

The Lagrangian of the axial CiDM model in [12] is given by 

^ = ^SM + ^* + ^Gauge + -^Higgs (1-1) 

where 



, ,. , - _ 1 p2 I l^pM»^D 

-Gauge — 2 -"-^ 2^ fJ-f 
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The fermion sector has an C/(2)ieft x C/(2)right chiral flavor symmetry, broken down to 
a vector-like U{1)h x U{1)l by Yukawa interactions with a dark Higgs boson, 0. The 
vacuum expectation value of 0, (0) = v^, causes the Abelian gauge field to acquire a 
mass tua^ — 2-\/2gfd^0 and the fermions to pair into a light and a heavy Dirac fermion. 



(1.4) 



with masses mi = yiv^ and rriH = ynVfi,, respectively. In this minimal model, rriH 
is near the electroweak scale while the other quark has a mass at or beneath the 
confinement scale, Aa ~ C(100 MeV - 10 GeV). 

In analogy to the Standard Model without weak interactions, the U{1)h x U{1)l 
flavor symmetry renders the lightest mesons and baryons charged under this symmetry 
stable. The lightest stable bound state is a meson containing a single and a. ^l, 
which will be denoted as tt^- This dark meson is the dark matter candidate in [12]. 
Because the constituents are fermions, TTd is paired with a vector state, p^. The spin- 
spin interactions of the constituents generate a hyperflne splitting 



hyperfine 
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where at is the SU{Nc) 't Hooft couphng, and the parameters are chosen such that 
the sphtting is at the scale (9(100 keV) suggested by DAMA/LIBRA. The hyperfine 
structure is described in more detail in Sec. 2.3. 

Another key feature of this model is that the dark gauge boson, A^^, couples axially 
to the dark quarks. The coupling between the A'^ and the dark mesons are constrained 
by parity and all leading order scattering channels are forbidden but the Hd — pd transi- 
tion [12]. A'^ mixes with the Standard Model photon and mediates dark meson/baryon 
scattering off SM nuclei. In particular, ttj up-scattering can explain the DAMA/LIBRA 
annual modulation signal. CiDM direct detection phenomenology is discussed in detail 
in [12,45,46]. 

The mass of the dark Higgs, </>, is radiatively unstable and introduces a second 
gauge hierarchy problem to the dark matter - Standard Model theory. The solution to 
the Standard Model's gauge hierarchy problem may also solve this new hierarchy prob- 
lem. Supersymmetric extensions of these models may solve both hierarchy problems 
at once, and may introduce new phenomena into the theory. For instance, if the only 
communication of supersymmetry breaking to the dark sector occurs through kinetic 
mixing, then the dark sector may be nearly supersymmetric, resulting in nearly super- 
symmetric bound states [47,48]. These susy bound states may have different scattering 
channels and the phenomenology may be different than minimal CiDM [49]. 

2. Dark Matter Synthesis Prelude: Hadron Spectroscopy 

The origin of the dark matter abundance in CiDM is non-thermal, possibly a baryogenesis- 
like process that generates a non-zero H — L number, see [66-77] for examples. Most 
often, particularly at large Nc-, no net dark baryon number is generated, resulting in 
a zero H + L number. This article addresses how the residual H and L asymmetry is 
configured at late times, i.e. whether the states are in dark mesons or in dark baryons. 

At temperatures above m^/, there is a thermal bath of gluons and free quarks and 
anti-quarks. When the temperature drops below m^^, the heavy dark quarks rapidly 
annihilate, and only the asymmetric abundance of H remains. The non-relativistic 
heavy quarks are much heavier than the confining scale and can form Coulombic bound 
states through antisymmetric color channels. Light L antiquarks have binding energies 
comparable to the confinement scale and therefore never bind at temperatures T ^ Ad. 
At confinement, the heavy quarks are dominantly screened by light antiquarks to form 
SU{N^) singlets. 

The evolution of the H and L asymmetry after confinement is dominantly deter- 
mined by the spectroscopy of the bound state systems. After the heavy quark mass, 
the largest energy scale in the problem is the Coulombic binding energy between the 
heavy quarks. The energetics of these heavy quark configurations guide the synthesis 
of baryon states through hadron binding reactions. The analysis of Sec. 3 illustrates 
the qualitative behavior of the synthesis of the H and L asymmetry in the early Uni- 
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verse. The fraction of dark mesons synthesizing into dark baryons depends sensitively 
on the masses of the hght hadrons and the results are parameterized by the mass of 
the lightest hadron. 

2.1 Coulombic Binding Energies 

The Coulombic binding energies of collections of heavy quarks is relevant for dark 
matter synthesis. Collections of uh < Nc heavy quarks can form quasi-Coulombic 
bound states that are deeply bound relative to the confinement scale when rriH 3> Aj. 
The deepest bound states of uh heavy quarks are antisymmetric color configurations. 
This section computes the binding energy of multiple H bound states using a variational 
approach. 

The Coulomb potential between two heavy quarks, i and j, in an antisymmetric 
color configuration is given by [50-53] 

V.,(r.,) = ^^C(2, 1, 1) = _^ (1±^, (2.1) 
where g is the ordinary gauge coupling and the running 't Hooft coupling is defined as 

"t(/i) = — ^ — (2.2) 

and C{ni,n2,n3) is the following combination of Casimirs for nj-rank antisymmetric 
tensors of SU{Nc) 

C{ni,n2,ns) = ^2(721) - ^2(712) - 6*2(^3). (2.3) 
The Casimirs for rank n antisymmetric tensors are 

C,(n) = G(n) = ^^<^^(l + i.) (2.4) 

where n = Nc — n in the above expression. 

Notice that the potential in Eq. 2.1 for combining fundamentals into antisymmetric 
tensors is attractive but Nc suppressed for fixed 't Hooft coupling in contrast to a quark- 
antiquark potential. 

The confinement scale, Aj, is determined by the 't Hooft coupling evaluated at a 
scale II = niH 

= ("(11 -2/^>.(m„))- P-^' 

for the minimal model with a single light flavor. 
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The total Hamiltonian for a system oi hh heavy quarks is 

2 1 

1=1 l^j 

In the ground state configuration, the color indices of the quarks are antisymmetric 
with all other quantum numbers being symmetrized. In particular, all the quarks will 
have the same spatial wave function in the ground state. The ansatz for the bound 
state system will be 

vl>("-)(xi,...,x„,) = nV^(^^), (2.7) 

i=l 

where ip{x) is a hydrogen-like ground state wave function of the form 

^(^) = (^) ' e-^-/^ (2.8) 

where /i,,, is the variational parameter. Evaluating the expectation value of 1-L gives an 
estimate of the binding energy 

(^(/i*)) = nn - ^^\^p inn - . (2.9) 

Minimizing ("H) with respect to gives the inverse Bohr radius for the hh heavy quark 
system: 

^i,{nH) ^ {UH - l)'^^mH{l + iV-i), (2.10) 
and binding energies of 

EB{nn)=nn{nn-lfEs, , = ' ^(1 + iV-^ (2.11) 

where the running of the 't Hooft coupling has been neglected in the differentiation. 
For the case of = 2 the binding energy reduces to 

precisely the analogue of the Rydberg constant for the potential in (2.1). 

The estimates of (2.11) are crucial for two dominant rearrangement processes that 
occur 

Vl>(") _^ -qjin') _^ v]>(n+n') _^ ^ ^(n) _^ ^(n') _^ ^{n+n'-m) _^ gy(m.)_ (2.13) 
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1 and 2 — 7- 2 processes. 



Ignoring the effects of running on the 't Hooft couphng, the energy released in these 
reactions is 

EB{n + n') - Esin) - Eein') - m^d 

Eb (n + n' -m) + Eb (m) ~ EBin) - Eb (n'). (2. 14) 

The energy release from the first few reactions are tabulated in Table 1. The first 
reaction releases the least energy and is a potential bottleneck in the chain reaction 
that takes N^^^^^ 

Eqs. 2.11 and 2.14 are the main results of this section. They will completely deter- 
mine whether heavy quark synthesis into bound states takes place before confinement. 
The role of light states (in particular light quarks) during the unconfined phase is irrel- 
evant due to the large gluon entropy at such temperatures. After confinement, when 
the ^("■^^ states are color-neutralized by light quarks, the estimates of (2.11) and (2.14) 
are still expected to hold, since they physically correspond to deeply bound states of 
heavy quarks and are typically an order of magnitude larger than Ad and the binding 
energies of light quarks. The energetics of the rearrangement reactions (2.13) after 
confinement will be completely determined by the results of (2.11) and the parametric 
dependence on the mass of the lightest state in the spectrum, that is typically present 
among the final states (more about that on Sec. 3.3). 

2.2 Hadronization of Bound States 

There are two ways of constructing color singlets from Uh heavy quarks in an antisym- 
metric color configuration: mesons and baryons. For mesons, the color is neutralized by 
Uh light antiquarks in an antisymmetric color configuration, while for baryons the color 
is neutralized by [Nc — uh) light quarks in a color antisymmetric state. In QCD, the 
analogues to these multi-heavy quark systems would be di-heavy baryon and tetraquark 
states of the form 

ccq,cbq,bbq ccqq' , cbqq , bbqq' (2.15) 

where q and q' are light-flavored quarks [54-59]. Only the ccq has been observed, but 
there is an ongoing program to discover the rest of these states. 
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In the dark sector of CiDM, the corresponding states will be denoted by ir^ " for 
meson states^ with rifj heavy quarks and by for baryon states with uh heavy 

quarks. For convenience, B^^"^ = Bh and 5^'^^ = B^. 

Due to the difference in the number of light quarks necessary to hadronize the color 
of vr^"'^'' and B^'^"\ these states will have different masses. Using a constituent picture 
of the light quarks, the masses of mesons and baryons are given by 

m (n„) = UHrriH - EsinH) + mesinn) + nispininn) 

mj^(„„) = UHrriH - EsinH) + mefr(A^c - Uh) + ms^in{Nc - uh). (2.16) 

In Eq. 2.16, mefr(?2L) parametrizes the constituent mass of the light quarks confined 
in a bag of size Ad 

mcairiL) = riiil + mo(nL))Ad. (2.17) 

In principle this leading order effective mass could be n^-dependent, due to changes in 
the light quark wavefunctions. For simplicity, this article will take mo('^L) = 0. 

The term mspininn) in Eq. 2.16 arises from the spin-spin interaction amongst the 
light quarks and a constituent model of this interaction gives 

mspi„(ni) oc C(2, 1, 1) Y,{S^ ■ S,). (2.18) 

Evaluating the expectation values above for a totally spin-symmetric state, one finds 

mspiniriL) = mi(nL)ni^^^^(l + iV,^')Ad. (2.19) 

Using the A;, — S;, system in QCD one infers: mf''(2) = 1.14. This article will use 
mi(?T,i) = 1.0. 

Both the constituent masses and spin-spin interactions cause the mesons to be 
lighter than the baryons when uh < and the baryons to be lighter than the 

mesons in the complimentary case. 

The most relevant reactions in the early Universe that synthesize stable B^""^ and 
vr^"^ states have rates that are exponentially sensitive to total binding energy differences. 
The estimates of the binding energies have an uncertainty of 0(Ad) that has been 
absorbed into the unknown constants in Eq. 2.16. In practice, the binding energy 
differences are larger than Ad by roughly an order of magnitude so that the Coulombic 
spectroscopy dominantly determines the synthesis of tt'^^ into other species of dark 
hadrons. 

^TT^"^ refers to all spin configurations of n and n dark quarks. Specifically, tt^^-* refers to 
both TTd and pd- 
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2.3 Hyperfine Structure Spectroscopy 

The stable dark hadrons have a large degeneracy of their ground states from the sup- 
pressed spin-spin interactions of the heavy quarks with the light quarks. Both in dark 
mesons and dark baryons, same-flavor quarks have antisymmetrized colors and sym- 
metrized spins. Therefore, heavy quarks are in a Sh = nH/2 spin conflguration while 
hght quarks are in a S'^ = nH/2 spin conflguration for mesons and Sl = {N^ — nH)/2 
spin conflguration for baryons. The resulting range of spins for dark mesons and baryons 
is, respectively. 



< J^^i^m) < riH 



< Jb(-h) < In,. (2.20) 



Hence, the degeneracy of a bound state containing uh heavy quarks (either mesonic or 
baryonic) is equal to nn + i- Spin-spin interactions break this degeneracy and introduce 
hyperflne splittings, resulting in the lowest spin conflguration being the lowest energy 
conflguration. Such splittings can be very suppressed with respect to the dark matter 
mass and offer a natural mechanism to generate the (9(100 keV) scale suggested by 
DAMA. 

The spin-spin coupling of heavy and light degrees of freedom splits the ground state 
degeneracy of hadrons. This is studied for Standard Model b hadrons in [60-62]. The 
quark chromomagnetic moment is suppressed by its constituent mass, m^,, and is given 
by 

ra»A, ^^^^^^ 




where the coupling is evaluated at the effective Bohr radius for the entire, color-neutral 
bound state. For m <^ Aj, m^, — )■ A^, at — )■ 1, and k is introduced to flx the relationship. 
The flrst order correction to the energy levels due to spin-spin coupling of heavy and 
light degrees of freedom is 

^Hyperfine ^ -C(0, , Tl/^ ) |^/^ | ■ S^l) |7/;(0) 1^, (2.22) 

where Sh and Sl are the spin operators for the collection of heavy and light quarks, 
respectively. The energy splittings can be evaluated using 

Sh ■ Sl = iji^"^ ^ ^li ^ ^li)-: (2.23) 

where S = Sh + Sl- The 5"!^ and Si terms do not induce mass splittings. The color 
factor C(0, uh, uh) in Eq- 2.22 is deflned in Eq. 2.3 and is given by 

C(0, riH. fin) = riniN, - uh) + (2-24) 
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Note that the color expression above apphes to both mesons and baryons. 

Finally, ijj{0) in Eq. 2.22 is the ground state wave function (i.e., n = Z = 0) for the 
light quarks evaluated at the origin and is roughly 



2^ 1 I l/a% = {atrriL)^ > Ad 



(2.25) 

47r [(^Ad)^ TTiL < Ad. 

For states with multiple light quarks, i.e. all states but vr^^'' and there could 

be significant multi-body effects that could change the wave functions. In particular, 
QCD baryons have a hyperfine splitting smaller than that of mesons by a factor of 
roughly three. For bottom hadrons 

C2 c2 / b . 2 

-f -f '-={^^] ^3.23. 2.26 

Remarkably similar ratios hold for charm hadrons and even strange and light flavored 
hadrons. This is evidence for hl dependence in k and can be interpreted as mo from 
Eq. 2.17 being 

rhoiuL) = — - 1. (2.27) 

This dependence in should be universal and applicable to both baryons and mesons. 
Combining the results for ^ Ad gives 

£^Hyperfine ^ < {Sh " S l) Uu f 1 - ^ (l + 0{N:')) (2.28) 

where rii = uh for vr^"'^'' states and ul = Nc — uh for B^"^"^ states. The energy splitting 
between the ground state and the first excited state for dark mesons is 

AE,„,)^«:2 n^/'l-^')^, mi«Ad. (2.29) 
For baryons the hyperfine splitting of the ground state is 



nn_ 



A 



+ 1]^, mi<Ad. (2.30) 



Fitting expression (2.30) to the hyperfine splitting of B and B* in the Standard Model, 
one finds 

kI ^ 1.35 (2.31) 

for Aqcd = 350 MeV and rrify = 4.5 GeV. Throughout this paper this constant will be 
taken to be ki = 1.0. 

QCD indicates k„ < ki, potentially making the smallest hyperfine splitting be 
in the vr^"''* meson, rather the vr^^'' meson. This opens up an alternate explanation of 
DAMA: that it is the multi-heavy quark mesons that are responsible for scattering. 
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Decay Mode 


Vd 
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4£ 
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Ai 




1— 


2£ 



Table 2: Examples of light mesons 



2.4 Light States 

Just like hybrid mcsonic and baryonic bound states, unstable states are in thermal 
equilibrium with the thermal bath at early times. The unstable spectrum consists of 
states that carry no net H and L quantum numbers. There are no Goldstone bosons, 
since this is a single flavor theory at confinement and the global U{1)f is anomalous. 
Therefore the mass of the lightest state is of the order of the confinement scale and could 
either be a light meson or a glueball state [63-65]. This section describes the unstable, 
0{Ad) mass states, estimating their lifetimes and briefiy describing their cosmology. 

Long hved states are of potential concern to Standard Model BBN. We will show 
here, however, that a general depletion mechanism causes the abundance of quasi-stable 
states to be negligible by the time BBN begins. Among the long lived states is the 
Vd ~ ^Ljb'^L meson, which decays through the operator 

A. = ~VdFd,.Fr (2.32) 

to four Standard Model leptons through two off-shell ^^'s- The resulting decay rate 
scales as [33]: 

. / 1 -ymt 



giving Tjd a typical lifetime of 



"Vd 



Scalar glueballs with J^'-' = mix with rjd mesons and inherit the same decay 
channels. J^^ = O'^^ scalar glueballs, on the other hand, mix with and can decay 
to either on-shell or off-shell A^s or light dark mesons. Their lifetime is expected to 
be no longer than the hfetime of the 0~+ states. 

Among the short-lived states are uj^ vector mesons and J^^ = 1 glueballs that 
can decay through mixing with oud- Their estimated lifetime is [33]: 



— e aao 

3 / m. 



r^d - ( -i^^'^^D ) (2.35) 
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The existence of short hved states in the spectrum is a leaky bottom mechanism 
for quasi-stable flavorless hadrons. These long lived states are depleted by scattering 
into shorter lived ones, e.g., rj^ + X ^ + X' , with cross sections set by A^: 

(^^),d+x^a.d+x' ^ 1/Ad- (2.37) 

The residual abundance of quasi-stable particles by the time Standard Model BBN 
begins is 

= -fm,^ ~ 10-'' GeV, (2.38) 

for m^^ ~ Ad. 

The most stringent constraints for a late decaying relic X with hadronic branching 
fraction _Bhad = 1 are set by '^He/D [79]: 

Cx = fmx < 10-'' GeV. (2.39) 

The relic abundance (2.38) is several orders of magnitude below the upper limit (2.39). 
This is a conservative bound because the decay products of the dark hadrons are 
predominantly leptons or photons, which alter the primordial abundances of elements 
less than hadronic decays. Therefore, light states have a negligible effect on BBN. 



3. Dark Matter Synthesis and Evolution 

The evolution of the H and L asymmetry in the early Universe is divided into three 
stages ordered by the temperature relative to the confinement scale: above, at, or be- 
low. Starting at temperatures close to mn, the thermal abundance of dark if-quarks 
is exponentially suppressed and only the non-thermal component is left to synthe- 
size into composite states. At temperatures above the confinement scale, perturbative 
techniques are applicable to estimate the cross sections for rearrangement reactions 
because rriH ^ Ad. For temperatures above confinement, the light quarks can be ig- 
nored throughout, and included during confinement to screen the color charge of the 
heavy quark bound states. Below confinement the light quarks can be treated as spec- 
tators. Assuming that approach is justified, the computation is reasonably accurate up 
to 0(l)-factors. 

3.1 Preconfinement Synthesis 

At high temperatures, T ^ Ad, the non-Abelian gauge dynamics is unconfined and 
quasi-perturbative. Light quarks are always relativistic and do not form bound states. 
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Below temperatures T ~ 0{mH), heavy quarks can form Coulombic bound states 
through antisymmetric color channels (as described in Sec. 2.1) by reactions of the 
kind 

^(n) + ^K) ^ ^(n+n') + ^ , Q = Ql^r} (^X^ = 0) (3.1) 

where \E'*^"-* is an antisymmetric bound state of n heavy quarks {2 < n < Nc), g is a, dark 
gluon and the energy releases, Q"^} and Q'^nn' given in Eq. 2.14 and tabulated in 
Table 1 with mx^ = rrig = 0. 

The capture cross section for reactions (3.1) is given by [80], 

47r^ Qt Q'nn/ 

where is the reduced mass of '^^"'^ + \Ef("'), (3/i„.„/T)^/^ accounts for Sommerfeld 
enhancement, and at(/i*) is the 't Hooft coupling evaluated at the Bohr radius of the 
bound state yU.* = atmn/^Nc. 

Each reaction of the type (3.1) contributes to the Boltzmann equation for the 
number density of n„, as^ 



— + SHUn = {avjr 



n„n„/ - n„+„/ ( exp(-Q„ „, /T) 

ZTT nin+n' J 



+ 



The exponential factor in the Boltzmann equation above accounts for the large gluon 
entropy that prevents heavy quark bound state formation down to temperatures 

- - ^^^2^ (3.4) 



In(FH) 



where Yu = nujs. Plugging in Eq. 2.12, we see that strong gluon entropy dissociation 
is effective down to temperatures when confinement takes place for: 

< —\\w.{Ju)\^d (3.5) 



GeV ~ «i(/x,)^ 



In 



GeV J 



21 



A. 



(3.6) 



For Ad > C(100 MeV) and A^'c = 4 that corresponds to m/^ < C(10 TeV). Thus 
for the parameter space favored by DAMA/LIBRA, there is no pre-confinement bound 
state formation. 



^with appropriate factors of 2 included whenever n 
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3.2 Synthesis at Confinement 

It was shown in Sec. 3.1 that dissociation due to relativistic dark gluons inhibits bound 
state formation down to temperatures T ~ Aj. At confinement, this dissociation shuts 
off since the theory acquires a mass gap. The light quarks and gluons form massive 
dark hadrons and proceed to decay promptly (Sec. 2.4). 

The formation of heavy quark bound states is suppressed since those are dilute at 
the time of confinement and there is a strong nucleation of light quark-antiquark pairs 
that screen the color charge of free heavy quarks. So, confinement will preferentially 
lead to formation of dark mesons over higher-n// dark hadrons. 

It is possible, however, to estimate the abundance of uh = 2 dark hadrons formed 
at confinement by simple combinatorics. If two heavy quarks are apart by a distance 
smaller that A^^ at confinement, the light-quark screening effect will not operate and 
the two heavy quarks will bind. Computing the probability that that will happen gives 
an estimate for the ratio of uh = 2 dark hadrons over tt^ dark mesons produced at 
confinement 



~ —IT ~ lU . 



(3.7) 



These bound states of two heavy quarks will either be screened by two light anti-quarks 
to form a vr^^^ dark meson, or by {Nc — 2) light quarks to form a i?*-^-* dark baryon. The 
later is Boltzmann-suppressed over the former for A'^c > 4 since it has {Nc — 4) more 
light constituents. 

The formation of dark baryons at confinement depends sensitively on Tconf- Using 
the leading term from Eq. 2.17 

^^expf-(Ar,-2n)-^y (3.8) 

n^ir,) \ I conf / 

Using QCD as a guide, Tconf — 0.3 Aj, leads to 

^ = 3 X 10-^ (3.9) 

for Nc = 4. The Uh = 1 dark baryons can annihilate with Bi into + X^- The 
residual B^^^ abundance that does not annihilate will quickly synthesize into Bh- 

3.3 Postconfinement Synthesis 

After the transition from the unconfined to the confined phase, the dark quarks hadronize 
dominantly into vr^^'' dark mesons with a suppressed fraction in fi^^^ dark baryons. 
Heavier dark hadron synthesis occurs through 



i'^ + ^i'^ ^ + 4°^ > i?W + 4^U5(2) + vrr, (3.10) 
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where tt^^^ is the hghtest dark hadron, such as a glueball or hght meson, with mass 

(2) 

miight — 0(Ad). Once vr^ mesons start forming, a chain of reactions can occur that 
will ultimately lead to formation of dark baryons Bh. 

As an illustration consider an A'c = 4 theory. The reaction chain down to heavy 
dark baryon synthesis is: 





+ 




Q 


—'^Ebo — ""^light 


(3.11a) 








Q 


=2Ebo + (mcff(2) - miight) 


(3.11b) 


-J" 






Q 


= 10Ebo — "flight 


(3.11c) 








Q 


=8Ebo 


(3.11d) 


-J" 




-^B„ + Bi 


Q 


=2AEbo 


(3. lie) 




+ 4^' 


-^B„ + Bi 


Q 


=32Ebo- 


(3.11f) 



As a combination of the mass gap in the theory and the fact that the heavy quarks in 
TT^ are not particularly deeply bound, the first reaction (3.11a) becomes endothermic 
for a large fraction of the parameter space. On the other hand, the B^^^ reaction 
(3.11b) is always exothermic; therefore, dark baryons formed at confinement process 
more efficiently than dark mesons. The results listed in Sec. 3.4 assume that all hybrid 
dark baryons efficiently process into B^- 

Reactions that have hybrid baryons in the final state 

ttJ") + vrf ) ^ + 5^ (3.12) 

are energetically disfavored compared to their mesonic counterparts, 

4") + vrf)->4"+"') + vr(°), (3.13) 

and have negligible contribution to the dark matter synthesis. 

Heavy quark binding in the rearrangement reactions (3.11) requires large momen- 
tum transfer Pmin = ruHV ~ \/2miiT ^ A^, and hence it is expected that the cross 
sections are controlled by perturbative heavy quark dynamics, as it was during pre- 
confinement (3.3). However, there are two important differences between pre- and 
post-confinement processing: 

1. There is no Sommerfeld enhancement because the lightest degrees of freedom are 
heavier than the typical momentum transfer in such reactions: mught > Pmin = 
mnv. This amounts to a reduction of the post-confinement cross section by a 
factor of (at/f) relative to the pre-confinement cross section. 

2. If miight is heavier than the vr^ binding energy, the first reaction (3.11a) is en- 
dothermic^ which further suppresses post-confinement synthesis by an additional 
Boltzmann factor e~*^™'"sht-£^B(2))/T ^.j^g thermally averaged cross sections. 
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Therefore, post-confinement cross sections for heavy hadron processing compares to 
pre- confinement heavy quark binding (3.3) as: 

Mpost=(^^^) 'e-(-'---^^-)/^(a^)p,e. (3.14) 

For the range of parameters where these reactions are endothermic, processing will 
be already frozen out by the time of confinement. This demonstrates that no heavy 
baryons are formed for mught ~ > Eb{2), or: 

A.^(^)'™». (3.15) 

However, that is not the only possibility. As one explores other ranges for the confine- 
ment scale Ad and the mass of the lightest state mnght, the full numerical solution to 
the Boltzmann equations reveals that the synthesis of CiDM can allow for much richer 
range of compositions detailed in Sec. 3.4. 

3.4 Synthesis Results 

Fig. 1 illustrates the CiDM synthesized spectrum as a function of the Ad — "mught 
parameter space for the case iV^ = 4. A general classification into five regions of 




Figure 1: Left: Dark matter synthesis as a function of the confining scale in the dark sector 
Ad and the mass of the lightest state mught for an A'^c = 4 dark sector. The inelastic splitting is 
6m = 95 keV. The different regions are described in the text and in Table 3. Right: Synthesis 
for Ad fixed at 190 MeV and 5m = 95 keV. 



different qualitative behavior is possible: 
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Region I: Complete Synthesis is efficient and the CiDM composition is dominated 
by heavy baryons Bh- DAMA's signal may arise from inelastic scattering of 
the highly suppressed dark meson components, particularly of tt\ if 
decreasing function of n^. The viability of this region requires that the heavy 
and light baryons are not visible in direct detection experiments. That will be 
further discussed in Sec. 5. 

Region II: Nearly Complete Heavy baryons are the dominant component, with a 
few percent of the dark matter density in the form of dark pions vr^ , n^^^ and irf^ . 
This region is not as extreme as Region I for CiDM. As with Region I, DAMA's 
signal may arise from scattering of the sub-dominant dark meson component. 

Region III: Incomplete Synthesis results in a democratic abundance with compa- 
rable mass densities for all for states ti^^\ t^^^\ v"^^^ and B^- Here the dominant 
number density of dark matter is the vr^^^ state and there are other components 
of the dark matter to discover. 

Region IV: Arrested The complete synthesized and unsynthesized components, Bh 
and TTd, share comparable mass densities, with a few percent in the form of exotic 
pions, TT^^^ and vr^ . The ffist step of the synthesis chain, tt^^V^^^ — )■ tt^^V^"'' is the 
bottleneck much like deuterium formation slows BBN in the Standard Model. It 
only occurs for a brief period, but once the ti\ has formed, it processes quickly 
into Bh- 

Region V: Inhibited The ffist step of the synthesis chain is strongly supressed and 
the CiDM composition is dominated by vr^^^ Region V is the cosmology taken 
in [12]. The heavy baryon component mostly arises through the primordial B^^^ 
formation described in Sec. 3.2. 

A quantitative description of the abundances in each one of these regions is summarized 
in Table 3. 



Region 


P„{i)/Pdm 


P„(2) / Pdm 


P_(3)/pdm 


PBh/Pdm 


I 


10-4-0.1% 


10-4 - 0.2% 


10-3 - 0.9% 


> 99% 


II 


0.1% -4% 


0.2% - 5% 


0.9% - 11% 


80% - 99% 


HI 


4% - 57% 


5% - 24% 


11% - 17% 


9% - 80% 


IV 


57% - 99% 


<5% 


<5% 


1% - 30% 


V 


> 99% 


< 10-5 


< 10-5 


< 1% 



Table 3: The relations on the fractional mass densities that define the regions of dark 
matter synthesis in Fig. 1. 
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The region of interest for DAMA/LIBRA realizes the full variety of dark mat- 
ter synthesis possibilities. In all cases, Hd is a good candidate for explaining the 
DAMA/LIBRA signal through inelastic kinematics, but signals from other dark matter 
components may be important. This will be discussed further in Sec. 5. 



Sec. 3 demonstrated that for a large range of the parameter space favored by DAMA/LIBRA 
[12,33,46] CiDM has a dark matter halo dominated by uh = 1 mesons, with subdom- 
inant components in other configurations. This section describes how the spin 1, pa 
dark mesons become depopulated. The discussion of this section is couched in "Region 
V" taking in the minimal CiDM model of [12] with a vTd dominated halo. This analysis 
can be generalized to any of the synthesis scenarios considered in Sec. 3.4 and the 
results do not significantly change. 

The nearly degenerate hh = 1 meson spin states and pd are equally populated 
at high temperatures. When kinetic decoupling between these two spin states takes 
place, the excited dark matter fraction, Up^/nTr^, is frozen in. If the excited state pd is 
cosmologically long lived, the fractional number density will be constrained from direct 
detection of pd — )■ tt^ de-excitation in nuclear scattering to be Up^/riT^^ < 10~^ [26,27]. 

These bounds are relevant whenever kinetic mixing of t/(l)d with hypercharge is 
the only channel for pd to decay. Then, the only kinematically allowed decays are to vTd 
plus photons or neutrinos. The decay to photons is highly suppressed (a loop-induced 
3-photon decay). The kinetic mixing decay to neutrinos is suppressed by an additional 
factor of eq^ /m\^ in the amplitude. The resulting pd lifetime is much longer than the 
age of the Universe in these cases [26,27]. 

4.1 Early Time Depopulation 

Down scattering constraints from long-lived excited states are a common challenge to all 
iDM models coupled to the Standard Model only through kinetic mixing of U{\)d with 
hypercharge. Kinetic decoupling of the spin states usually occurs before T ^ 1 MeV, 
leading to an 0{1) fractional number density np^/riT,^. Strongly coupled composite 
dark matter models exhibit an elegant way to depopulate the excited states. 

The reactions that keep the dark meson spin states in kinetic equilibrium have 
large cross sections set by the size of the composite 



where fo is an 0{1) constant that parameterizes the quasi-elastic scattering cross sec- 
tion. These reactions will freeze out when 



4. The Excited Fraction of CiDM 




d 



(4.1) 




~ 3e 



5m./T, 



(4.2) 



PdPd-S>'rdi"d 
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where 6m ~ (9(100 keV) is the energy sphtting between and vTd, is the freeze 
out temperature and H{T^) is the Hubble expansion rate at freeze-out. With ~ 
niHSm/ nf, Eq. 4.2 imphes that the primordial up-scattered fraction np^/n^^ ~ 3e~^"^^'^* 
is completely determined by the dark pion mass: 

^-Z"- -5xl0-.?f,^ (4.3) 



\\nnpjn^j ' VlOO GeV 

In other synthesis regions where the vr^^^ abundance is suppressed, e.g. Region 
I or II, the depopulation if pd is more effective because the dominant de-excitation 
interaction is 

Bh + Pd^ Bh + TTd. (4.4) 

This interaction has roughly the same cross section as (4.1), but Bh acts as a catalyst 
for de-excitation resulting in 

„.,.X-<.(r, . r..2eV(^)'^(^)(^) (4.) 

which severely depletes the pd population. 
4.2 Late Time Up Scattering 

After structure formation begins, the velocity of the dark matter increases allowing for 
late-time up scattering of ir^ into pd in dark matter halos. The excited state can then 
be repopulated to observable levels today. The up scattering reaction 

TTd + TTd -> Pd + Pd (4.6) 

is endothermic, so its cross section scales as 



1 6m 1 

The number of up scattering collisions since structure started forming, about Tstruct ~ 
10^'^ yr ago, is roughly: 

, -I / 100 GeV\ VlOO keV\ , 

with i^viriai — 10^'^. For mTTj ~ 100 GeV, a fraction 0(10""^) of the vTd have up-scattered 
during the entire age of the Universe. Up-scattering from reactions of TTd off nuclei 
through y4j-exchange can generate a similar pd abundance. 

Eq. 4.8 demonstrates that, even though up-scattering during structure formation 
does not surpasses the primordial excited fraction for m^^ ^ (9(100 GeV), it is the 
primary up-scattering effect for m,rd < (9(100 GeV), reaching the current upper bound 
np^/riT,^ ~ (9(10^^) for dark pion masses as low as (9(10 GeV). 
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4.3 Down Scattering Discovery 

Although the predicted CiDM excited fraction np^jn-,,^ is safely below the current 
upper bound, it is in the reach of discovery of the next generation of direct detection 
experiments. In order to quantify that, in what follows we compute the detection rate 
of the pd component as it down-scatters off of target detector nuclei. 

The down-scattering interaction is determined by the operator {2gAl ^A)T^AdvpApF^^-, 
where mixes with Standard Model hypercharge B^^ (see Eq. 1.2). Its differential 
rate is 1/3 of the differential rate for up-scattering [12], 

Here Eji is the nucleus recoil energy, mjq and Z are its mass and atomic number, 
respectively, and frci is the pd-nucleus relative velocity. The Helm nuclear form factor 
\F{E^\^ accounts for loss of coherence scattering off of the entire nucleus at large 
recoil, and it is given by 

i^(^^)|2^^M^y^_,.i,.^ (4.10) 

where s = 1 fm, Tq = y/r"^ — Ss^, r = 1.2y4^/^ fm, and \q\ = y/2mj\iEji is the momentum 
transfer. Finally, f^s is defined by: 



(4.11) 



/eff ^A, 

where niA^ is the mass of the f/(l)d gauge boson. 

The differential rate per unit detector mass is given by: 

dR n„, p° f ^ . da . 

dF- = ^^^rt- / ^ ^ ^4-12) 

n^^ rn^^rriN Jy^^^.^ dh^ 

where is on the order of the local dark matter density, p^y^ ^ 0.3 GeV/cm^, since 
we are assuming here that vTd is the dominant halo component. The Standard Halo 
Model (SHM) velocity distribution function f{v,VE) with a functional form 

f{v, ve) = U (^e-^^^+'^^^'/^'o - e-^'-/^o) e(t;esc -\v + ve\) (4.13) 

will be used in this article, with the following benchmark halo parameters: local escape 
velocity fesc = 550km/s, local velocity dispersion vq = 270km/s, and the velocity of the 
Earth in the galactic frame as in [25]. Moreover, the minimum relative velocity for 
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Figure 2: Predicted number of events per kg-yr of exposure due to pd — )• vra down-scattering 
off of a Ge target, and detector sensitivity in the 10-100 keV range. The gray region is 
excluded by the CDMS 2009 Ge analyzed data [83] at 95% C.L. 



Pd to inelastically down-scatter to vTd causing the target nucleus to recoil with energy 
Eji is given by 



^min(-E'ij) 



-5m + 



(4.14) 



where /i is the DM-nucleus reduced mass. 

Figure 2 displays the predicted rate at a detector with Ge as the nuclear target, 
with sensitivity in the range 10 — 100 keV, such as CDMS. The region favored by 
DAMA/LIBRA evades all current bounds, but the next generation of direct detection 
experiments, such as XENONIOO, should have sensitivity to discover it. The distin- 
guishing features of such down-scattering signals that can help disentangle it from other 
components are: (i) lower recoil energy threshold below which the signal vanishes, since 
the hyperfine splitting energy is converted to kinetic energy after de-excitation, and (ii) 
for the region where the signal peaks, around Er = 6m^^, there is a phase inversion 
in the annual modulation relative to the up-scattering rate {i.e., the pd T^d rate peaks 
on Dec. 2°*^) [82]. Nevertheless, the modulation fraction is of order of a few percent, 
as for elastic scattering, since most of the halo is kinematically allowed to scatter. 
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5. Baryon Elastic Component and Direct Detection 



The calculations of Sec. 3 indicate that it is possible to have significant synthesis of 
the TTd into baryons. When this occurs, there are equal numbers of heavy baryons and 
light antibaryons. This section considers the direct detection implications of baryons 
that can arise through synthesis in the early Universe. These baryons have large spin 
since they consist of single fiavor quarks, and the axial constituent coupling will lead 
to a small but potentially detectable elastic scattering channel. 

5.1 Heavy Baryon Scattering 

As discussed in Sec. 2.3, a Bh heavy baryon has no hyperfine splitting since all its 
constituents have the same fiavor. Hence its dominant scattering channel is elastic. If 
its constituents were vectorially coupled to the dark group t/(l)d, dark heavy baryons 
would have an upper limit on their density fraction set by null searches (such as CDMS 
or XENONIO) to be in the range 10~^ — lO"*^ (assuming that the inelastic component 
fits the DAMA signal). That would imply that CiDM would have to live in Region 
V of the Ad — miight parameter space as described in Sec. 3.3; all other regions being 
excluded. 

The picture changes when the dark quarks have a purely axial coupling to U{l)d- 
We shall demonstrate that in that case the U (l)d gauge boson couples to the spin of the 
heavy baryon in the non-relativistic limit. This leads to a suppressed scattering rate 
relative to the case where the baryon has a vector-like U{l)d charge, opening up the 
parameter space to allow for all processing regions described in Sec. 3.3. One remark- 
able consequence of this framework is the possibility that the dominant component of 
the halo is in the form of dark heavy baryons and what the DAMA/LIBRA experi- 
ment is detecting is a subdominant inelastic component of the halo. In this section 
we will make these statements quantitative and consider their consequences for direct 
detection. 

The fermionic axial current is given by 



where the index n sums over constituent fermions and refers to their axial charge. 
Taking the non-relativistic limit and integrating out the small components, the expres- 
sion above for Nc heavy quarks reduces to: 





n 



Bh ~ 




where Xb is the Ac/2-spin wave function of the heavy baryon bound state. 



XBh = ^1 ® ^2 ® ... (S)^Pn, 



(5.3) 



i.e. the direct product of the two component constituent wave functions ipn- The heavy 
baryon spin operator S^^ is given in terms of the constituent spin operators as 

^Bh = ® I2 ® ••• ® liv, + • • • + li ® ••• ® l7Ve-i ® (5.4) 
Note that its Casimir invariant is given by 

SL ^ ^ (^ + iV (5.5) 



2 V 2 

and it satisfies the following algebraic relations 

[SW' = '^'''^k {Sb,, = Inc{N, + ^W'- (5.6) 

The differential cross section for a dark heavy baryon scattering off a target nucleus 
can be computed using Eq. 5.2 and the kinetic mixing between f/(l)d and t/(l)y, 
recalling that this results in the proton picking up an effective ee charge under U{l)d- 
The resulting cross section is 

da N,{N, + 2)aZ' 1 ^ ,,2 



dE^ 12(iV, + 1) vi, 
and the differential rate per unit detector mass is given by: 



En\F{En)\\ (5.7) 



- F,,^^ [ d^v vmv,)^. (5.^ 



dER mBHiTiN dE 

Here, Fp^ = Pbh/Pbm fractional dark matter density in the form of heavy baryons 

and mBjf = Ncmu is the dark heavy baryon mass. The velocity distribution function we 
use is the same as in the previous section, v^am denotes the minimum relative velocity 
for a dark heavy baryon to elastically scatter off a nucleus with recoil energy Er 



V2mNEn , . 

VminiEB) = — , (5.9) 

where /i is the ms^-nucleus reduced mass. 

Fig. 3 illustrates how the number of events expected at a Ge detector with sen- 
sitivity in the 10 — 100 keVnr nuclear recoil range (such as CDMS) depends on the 
fractional heavy baryon density for tuh = 72 GeV, assuming that the inelastic compo- 
nent fits the DAMA/LIBRA signal [45]. Note that if the halo is highly dominated by 
heavy baryons, with only a fraction of 10~^ — 10~^ in other components, then detection 
of the heavy baryon component is around the corner [83]. 
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1000 




Figure 3: Expected number of events per kg-yr from elastic heavy baryon scattering off Ge 
for niH = 72 GeV and nuclear recoil range 10-100 keVnr, taking the inelastic halo component 
to fit the DAMA signal. Fp^ denotes the fractional dark matter density in the form of heavy 
baryons. The green region is excluded by the CDMS 2009 Ge analyzed data [83] at 95% C.L. 
The dashed line denotes the projected sensitivity of SuperCDMS for 30 kg-yr of exposure 
at Soudan mine, the dot-dashed line 300 kg-yr at Snolab and the dotted line 3 ton-yr at 
DUSEL/GEODM. 



5.2 Light Baryon Scattering 

Recently, CoGeNT has seen anomalous low energy events consistent with light dark 
matter (m ~ 7 GeV) elastically scattering [88-91]. The rate is slightly larger than 
DAMA's rate if channelling is included [28,92-96]. A candidate for such light dark 
matter inside CiDM are the light baryons that arise while synthesizing the dark pions in 
heavy baryons. There are two issues with this interpretation, the first is that CoGeNT's 
rate would oversaturate the light baryon rate at DAMA's detector obviating the need 
for inelastic dark matter. However, channeling at low nuclear recoil energies is uncertain 
and may not be as effective as naive extrapolations from higher nuclear recoil energies 
indicate [97-101]. If channeling is less effective, then DAMA's signal may still arise 
from inelastic scattering, leaving CoGeNT's signal unexplained. 

The second, more serious problem with explaining CoGeNT's anomaly with light 
baryons is that there are equal numbers of heavy baryons produced. These heavy 
baryons have approximately the same per nucleon cross section as the light baryons. 
Heavier mass dark matter candidates have tight limits on their elastic scattering. 

There are two possibilities to evade the heavy baryon limits. Heavy baryons can 
decay into a light baryons and other colored particles. These types of decays generi- 
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cally result from the interactions that generate the asymmetry responsible for the ir^ 
abundance. It is possible that the heavy baryons decay while leaving the vTd stable 
resulting in a dark matter sector dominated by dark mesons and light baryons. For 
instance, if the dominant operator that violates H — L number has charge 4, e.g. 

>C_ff_Lviolation = 'Jp^^^^^^y (5.10) 

then TTd will be stable, all other heavy dark hadrons will be unstable. Specifically, the 
baryons will chain decay as 

fM o\ \ Bt + X Nr even 

Bh 5(^=^2) + X ^ > < (1) ■ 5.11 

[4 +^L + ^ AT, odd ^ ^ 

While these decays are model dependent, it might be possible to correlate cosmic ray 
signals from dark meson decay or oscillation with CoGeNT's signal. 

Another method for evading the heavy baryon constraints is to have a non-zero 
H + L number generated resulting in an excess of B^s. Generating a significant dark 
baryon asymmetry for large Nc is challenging because the dark baryon operators have 
high scaling dimension and therefore are highly irrelevant. For = 3,4,5 it may 
be possible to generate a sizeable asymmetry, particularly in supersymmetric theories 
because dark squarks have lower scaling dimension than their fermionic partners. 



6. Conclusion 



In Composite inelastic Dark Matter, the (9(100 keV) mass splitting suggested by 
DAMA/LIBRA arises from dynamics inside composite states of a strongly interact- 
ing sector. Dark hadrons carrying non-zero flavor quantum numbers are stable and 
hence potential dark matter candidates. In particular, heavy flavor mesons offer a 
good candidate for implementing iDM. The thermal relic abundance of dark hadrons 
is too small to account for all the dark matter due to strong self-interactions, meaning 
that the dark matter abundance in these theories must originate from a dark flavor 
(or dark baryon) asymmetry. One of the purposes of this work was to investigate the 
cosmological evolution of the flavor asymmetry and how it determines the mesonic and 
baryonic abundances at late times. This question was addressed in the context of the 
minimal model of Sec. 1.1, where the iDM candidate is a spin-0 dark meson TTd that 
inelastically scatters to a nearly degenerate vector state pd- 

In a large part of the parameter space, TTd mesons dominate the abundance of dark 
matter. In other regions, dark baryons dominate the abundance; however, the residual 
TTd component interacts sufficiently strongly to give rise to a viable iDM scenario. In 
all cases, exotic dark matter components arise with novel elastic scattering properties 
- nuclear recoil events are suppressed at low-energy by a dark matter form factor. The 
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relative abundance of pd to is typically ~ 10~^, suggesting that pd down 
scattering off nuclei is a discoverable signal in the near future. Finally, we studied a 
variety of long lived meson and baryons states and found that BBN constraints on their 
decays are not severe. 
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